Draft version July 11, 2008 

Preprint typeset using WT^^ style cmulatcapj v. 10/09/06 



00 

O 
O 
(N 



MODELING KICKS FROM THE MERGER OF GENERIC BLACK-HOLE BINARIES 

John G. Baker, William D. Boggs \ Joan Centrella, Bernard J. Kelly, Sean T. McWilliams \ M. Coleman 

Miller ^, James R. van Meter 

Laboratory for Gravitational Astrophysics, NASA Goddard Space Flight Center, Greenbelt, Maryland 20771 

Draft version July 11, 2008 

ABSTRACT 

Recent numerical relativistic results demonstrate that the merger of comparable-mass spinning 
black holes has a maximum "recoil kick" of up to ~ 4000 km s~^. However the scaling of these recoil 
velocities with mass ratio is poorly understood. We present new runs showing that the maximum 
possible kick perpendicular to the orbital plane does not scale as ~ rj^ (where 77 is the symmetric mass 
ratio), as previously proposed, but is more consistent with ~ rf , at least for systems with low orbital 
precession. We discuss the effect of this dependence on galactic ejection scenarios and retention of 
intermediate-mass black holes in globular clusters. 

Subject headings: black hole physics - galaxies: nuclei - gravitational waves — relativity 
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1. INTRODUCTION 

Recently, numerical exploration of the radiative recoil 
"kick" of merging black holes has progressed consider- 
ably. In particular, efforts in this regard have led to sug- 
gested phenomenological formulae for the kick, largely 
based on post-N e wtoni an (FN) predictions such as that 
given bv iKidde^ (|1995D . w hich have proved surp risingly 
successful. For example, iGonzalez et al.l (|2007[ ) found 
that in cases of unequal masses {q = toi/to2 < 1) and 
no spin, a simp le^ modifi cation of the FN formula orig- 
inally found by iFitchettI (jl983l ) fits the numerical data 
quite well. For cases of spins perpendicular to the orbital 
plane (i.e. parallel wit h the orbital angula r momentum), 
a formula proposed byl Baker et al.l (|2007f) is also consis- 
tent with numerical data. This formula is loosely based 
on FN calculations, with spins perpendicular to the or- 
bital plane producing kicks in the orbit al plane. For spins 
with c omponents in the orbital plane, ICampanelli et aTl 
(|2007t l have proposed a formula, again derived from FN 
calculations, that agrees well with numerical results for 
equal masses. 

This last type of kick, which is perpendicular to the 
orbital plane, is of particular interest because its com- 
puted magnitude can be very large (up to thousands 
of kilomete rs per second). In th e current literature 
(specifically ICampanelli et all l2007l ) , the mass-ratio de- 
pendence is drawn from the leading-order FN approxi- 
mation. It is unclear whether this approximation is suf- 
ficient to predict the strong-field dynamics that presum- 
ably determines the kick. Indeed, hints of a deviation 
from this f orm are evident for mass r atio g = 1/2 in 
the runs of iLousto fc Zlochowed ()2007l ). Therefore, al- 
though the angular dependence of the proposed formula 
is consistent with symmetry argumen ts, which are inde- 
pendent of the s trong-field dynamics (|Bovle et al.l 120081 : 
iBoyle fc Kesde n 2007), the mass ratio dependence of this 
formula is currently not well justified. 

Characterization of the dominant kick for unequal 
masses is especially important because, although the 
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largest possible kicks would eject the remnant from any 
galaxy, for astrophysical applications the distrib u tion of 
kick s peeds matters most. For example, iBonning et al l 
(|2007f ) find no evidence for quasars ejected from their 
hosts (although iKomossa et al.l |2008| may have seen a 
2650 km s~^ kick). If quasar activity is commonly in- 
duced by major galaxy mergers that lead to coalescence 
of supermassive black holes, the implications of this 
therefore depend in part on how frequently one expects a 
merger to allow ejection. Even for very large recoils, more 
than half of galaxies would still retain black holes at their 
cores (Schnittman 2007). However, the kick speed dis- 
tribution has a major impact on the hierarchic al growth 
of ma ssive black holes at redshifts 2; > 5 (e.g.. IVolonteril 
I2OOI . 

Here we investigate how the out-of-plane kick depends 
on the mass ratio, and find that, for mass ratios in the 
range 5 = 1 to g = 1/3 and spins Si < 0.2rrif, the kick 
drops off mo re rapidly with de c reasin g mass ratio than 
proposed by ICampanelli et al.l ()2007l ). Specifically, we 
find that a large body of numerical data on kicks are 
well represented by 
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where Vm is the mass asymmetry contribution, v± and 
U|| are the spin contributions that yield kicks perpen- 
dicular and parallel to the orbital angular momentum. 

r/ = q/ (l + q)'^ is the symmetric mass ratio. a|' is the pro- 
jection of the dimensionless spin vector di = Si /ml of 
black hole i along the orbital angular momentum, while 
and (t)i are the magnitude and angle with respect to 
some reference angle in the orbital plane of its projec- 
tion, cS^, into the orbital plane. $1 and $2 are constants 
for a given mass ratio. Here, A = 1.35 x lO^kms"^, 
B = -1.48, H = 7540 ± leOkms"!, ^ = 215° ± 5°, 
and K = 2.4 ± 0.4 x 10'''kms"^ This formula, simi- 
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TABLE 1 

Initial Configuration and Final Kick for Each 
Simulation. </>i(2) is the angle made by the spin 

vector of hole 1(2) WITH THE VELOCITY VECTOR 
OF HOLE 1, AS SHOWN IN FiG. [T] NUMERICAL RESULTS 

FOR THE KICK COMPONENTS Vm (WHERE AVAILABLE) 
AND D|| ARE SHOWN. KiCKS FOR EQUIVALENT SPINLESS 
RUNS ARE IN PARENTHESES. 
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lar in form to th at of rCamp ancUi ct al. (2007), synthe- 
sizes results f rom [ Gonzalez et al.. (,2007. ) for ([2]) and from 
iBaker et al.l ([20071) for (P and ^ ^ For $ and H we 
have f it available n umerical d ata from [H errman n et al.l 
([20071 ): iKoppitz et a l. (200'^: iBaker et al . (2007). The 
qualitatively new part, the factor of 7 7'^ in (^, replaces 
the fa ctor of 77^ originally proposed by ICampanelli et all 
([2007t l. and is motivated by new numerical evolutions 
presented here. 

2. INITIAL DATA AND METHODOLOGY 

We simulated the inspiral and merger of a range of 
spinning black-hole binaries, with mass ratios in the 
range 1/1.1 > q > 1/3. The initial configuration of 
momenta and spins is illustrated in Fig. [T] The parame- 
ters used in the numerical evolutions are presented in the 
first three columns of Table [TJ For these evolutions, the 
smaller hole (toi) has a dimensionless spin |ai| = 0.2, 
while the larger hole's spin is \d2\ = q^ldil. Both spins 
initially lie in the orbital plane, at angles 0i and (f>2 to 
the initial velocity of hole 1 (see Fig. [T]). 

To perform our simulations, we er nployed the Hahn- 
DOL e volution code, as described by IBaker et al.l ([20071 
[20083) • We chose initial coordinate separations of 7.0M 
for the q > 1/2 cases and 8.0M for the q — 1/3 
cases (where M is the total mass of the system) to 
yield between one and four orbits prior to merger; in- 
formed by PN theory ([Damour "etaI][20M ). we chose 
the corresponding momenta to minimize initial eccen- 
tricity. The finest grid spacing in all the runs presented 
here was hf = 3M/160. We also performed a single 
high-resolution simulation of hf = M/64 for the q = 1/2 




Fig. 1. — Configuration of black holes for all new simulations. 
The two holes' spins «(i,2) li^ initially in the orbital plane, at angles 
(f)l and {j>2 to VI , the smaller hole's initial velocity. 

case and found the kicks and all other relevant quantities 
agreed with the corresponding hf = 3M/160 simulation 
to within ~ 1%; additionally we performed a set of lower- 
resolution hf = 3Af/128 simulations in the q = 1/3 case, 
demonstrating consistency of the amplitude of U|| with 
the hf = 3M/160 simulations to within 6%. 

3. RESULTS AND DISCUSSION 

The recoil kicks resulting from the new simulations are 
in the rightmost columns of Table [1] Note that the in- 
plane kick agrees well with what would have been ex- 
pected given no spin («,„). This indicates negligible or- 
bital precession, which we also verified from the trajec- 
tories of the black hole centers. 

To conceive of plausible candidates for the mass scaling 
of I'll , we begin with the spin expansion and symmetry ar- 
guments of iBovle et all ([20081 ): [Boyle & Kesded ([2()07| ). 
For the spin configurations considered here, 

= Diq)ai cos (0i - ^q))~D{l/q)ai cos (02 ~ $(l/g)) 

(5) 

where D and $ are some functions of mass ratio q, and 
we note that $ must also depend on the initial sepa- 
ration. Further restricting ourselves to forms relatable 
to the factor of Si /mi — 82/ m2 app earing in PN cal- 
culati o ns of the kick, which inform ed ICampanelli et all 
([2007t l: lLousto fc Zlochoweil ([20071 ) and has been numer- 
ically well-verified in the equal-mass case, we substitute 
D{q) = qC{r])/{l + q) to obtain: 



Civ) 

(1 + 9) 



[qa^ cos((/)i - $1) - cos((/)2 - $2)] 



^ Note that in IBaker et al.l l|2007l ) , we used a simpler form for 
the zero-spin contribution, equivalent to ^ with B = 0. 



(6) 

where $1 = $(g) and <I>2 = '^(l/?)- Eq. [4| arises from 
the choice C(7/) = Krj^, with K = 2A x 10^ 

There are several ot her possibilities for the form of 
C{r]) in the literature. ICampanelli et al.l ([2007f ) assume 
C(?7) = 6.0 X lO'^rj^. In this cas e, ((/>i - ^1) and (^2 - ^2 ) 
are related to (6 — 9o) from ICampanelli et al.l ([2007f l. 
as well as the angle between the spin vectors implicit in 

\^2 " qdj^l- 

Another possibility arises from the known relation be- 
tween v\\ and the difference between the energy radiated 
in the ( l,m) = (2,2) and (2, —2) harmonics of the ra - 
diation ([Schnittman et all [20081 : iBriigmann et al.ll2007D . 
With no spin, these quantities are equal. With spin. 
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TABLE 2 

Maximum percent error resulting from various models 
of the kick, as distinguished by overall mass-ratio 
dependence. see equation (8). 



g I 

1/1. il 

1/1.3 
1/1.5 
1/2 
1/3 



K(a2ri^ +a4r;*) 



0.22 


0.23 


0.20 


0.75 


0.80 


0.78 


1.26 


1.31 


1.28 


15.57 


2.46 


1.41 


39.58 


10.98 


9.11 



we expect that ~ E22(peak)F ^ where E-^ 



^22(peak) 



is the 



peak power radiated in the (2, 2) harmonic, and F repre- 
sents the spin-dependent asymmetry between E22(peak) 

and i?2-2(peafe), i-e. F I - i?2-2(peafc)/-^22(pea*;)- For 

black holes with no spin, we have found that -E22(peafe) = 
a2rf -\- a4r]'^, where 02 = 0.0044 an d 04 = 0.0543, give s 
a good fit to the numerical data ([Baker et al.l l2008"bh . 
We do not expect spins orthogonal to the orbital angular 
momentum to change the scaling of the radiated energy 
significantly. If we further assume that the asymmetry 
factor F is independent of rj, which finds some support in 
PN analysis since to leading order P\\/E is independent 
of ry, then we hypothesize that C{rj) oc {a2rf' + a^rj'^). 

In Table [2] we summarize the agreement of various kick 
formulas with the numerical data. For each formula, 
which has the form of Eq. ([6]), we found the best <&! 
and $2, per mass ratio, according to a least-squares fit 
to the data given in Table [1] For each mass ratio, the re- 
sulting percent error is given for each model, maximized 
across initial angle. Referring to Eq. ([6]), the column 
headings Kij"^, K{a2ri^ + a/^-q'^) and Krf of Table [2] rep- 
resent choices for C{vi) that were tested, where in each 
case K has been chose n so as to reproduce the value of 
the formula of Campan elli et al.l (2007) in the equal-mass 
case. 

Now we consider th e agreement of our data with the 
C(ri) = Krf scaling of lCampanelli et al.l (j2007l ) (first col- 
umn of Tabic [2]) . The error of the best fit grows signifi- 
cantly with mass ratio, hence the mass-ratio-dependence 
of this formula is inaccurate. One might suppose that 
precession of the spins into the orbital plane could ac- 
count for this. However, the column in Table[T]shows 
that the in-plane kicks are close to those measured with- 
out spins (given in parentheses); hence this does not ex- 
plain the discrepancy in v\\ from the scaling. We have 
experimented with other values of K to resolve the dis- 
crepancy. For example, the maximum error of the 
model can be reduced to less than 10% for the q — l/i 
case, but not without increasing the maximum error of 
other mass ratios closer to unity to significantly greater 
than 10%. 

Since or iginal submission of this pap er, new data pre- 
sented by iLousto k, Zlochoweil (|2008D seem to indicate 
rj^ scaling, although the cases analyzed are complicated 
by considerable orbital precession. For example, their in- 
plane kicks are apparently at odds with previous formu- 
lae. It is possible that strongly precessing orbits require 



different fitting formulae, but this has yet to be settled. 

The choice C{vi) — K{a2rf -\- aiVj^), motivated above, 
fits the data much more successfully (second column 
of Table [2|) . Other scalings ca n be motivated throug h 
post- Newtonian-based analysis (Schn ittman et ani2008f l. 
However, a better empirical model was found to be 
C{rf) — Krf (third column of Table [2]). For now we 
consider this our best fit, and leave open the interest- 
ing question of how to accurately relate this prefactor 
directly to i?22- 

Our results affect the distribution of kick speeds, 
given various assumptions about the spin parameters, 
spin orientations, and mass ratios involved in coales- 
cences. This has particular application to the retention 
of the products of mergers of massive black ho les in 
the current universe (e.g., iBonning et al.l I2007D and 
electromagnetic signatu res of kicks (e.g.. [Shields et aTl 
l2007t iLippai et al.ll200^ . as well as coalescences in the 
early structure formation phase of reds hift 2: ^ 5 — 30 
(Merritt et al.' 2004; Bovlan-Kolchin et al. 2004; Haiman 
2004j : .Madau fc Quataert .2004; Yoo fc Miralda-Escudi 
20041 IVolonteri fc Pernal (20q5l: iLibeskind et all l2006f 



Micic et al.l 120061: IVolonteril |2007D T and for current- 
day mergers of intermediate-mass black holes 
(IMBHs), which might exist in dense stellar clu s- 
ters ( Taniguchi et all l200d. j Miller fc HamiltonI l2002blla 
Mour i fc Taniguchil boOM^ iMiller fc ColbertI I2004 



iGiiltekin et al .1 l2004Ll200a lO'Learv et al.l I2006L [2007 
Note that g = 1 to g = 1/3 is in the range of ra- 
tios expected for m ajor mergers of galaxies, and as 
ISesana et al.l ([2004D show, this range is expected to 
account for most massive black hole mergers in the early 
z > 10 phase of black hole assembly. 

Our new formula implies an important revision in our 
understanding of how easily IMBHs with M ^ 10^ — 
10^ Mq are retained in globular cluster s. A rich cluster 
has an escape s peed fesc ~ 50kms~^ ([Webbink| [l985). 
IGiiltekin et al.l ([2006D showed that the Newtonian kicks 
involved in binary-single interactions are insufficient to 
reach this speed if the IMBH is at least ~ 15 — 20 times 
more ma ssive than the obj ects with which it interacts. 
Using the lCampanelli et al.l (2007) formula, however, the 
maximum kick from gravitational radiation is Vmax = 
6x10'' kms^^?]^, implying that even IMBHs 30—35 times 
more massive than the black holes wit h which they merge 
could get ejected. Hollev-Bockelma nn et al.l (2007), fo- 
cusing on cases in which stars lose little mass through 
their evolution and thus can leave behind stellar-mass 
black holes with masses > 60 — 100 M©, use this to ar- 
gue that most IMBHs of even 1000 M© will be ejected 
from globulars. If instead stellar-mass black holes have 
masses ^ 10 Af©, a mass of at least 400 Mq would still 
be required to guarantee retention. 

In contrast, our new formula suggests a maximum kick 

of Wmax = 2.4 X lO^kmS"^?]^. Thus if 77 < 0.06, Umax < 

50kms^^. Therefore, an IMBH interacting with 10 M© 
black holes will stay in a rich globular if its initial mass 
is M > 170 M©, comparable to what is necessary for 
retention against Newtonian three-body kicks. 

Our results also have implications for whether merged 
supermassive black holes stay in their host galaxies. 
The figure of merit is the fraction of kicks that ex- 
ceed typical escape speeds from galactic centers (rang- 
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ing from roughly 500 kms^^ for a small spiral to 
2000 kms~^ for a giant elliptical), given assumptions 
about the distribution of spins and orbital orientations. 
The calculation of recor d for t his purpose is that by 
ISchnittman fc Buonanno (2007), who used a kick for- 
mula based o n effective one-body anal ysis and is different 
from that of ICampanelli et al.l (|2007f) : this formula un- 
derestimates the highest kicks. Table [3] compares the 
fractio n of kicks above 500 kms~^ and 1000 kms~^ us- 
ing the ISchnittman fc B uonanno (2007) formula (an un- 
derestimate). the ICampanelli e t al. (2007) formula (an 
overestimate), and our r esults . It is clear that the 
ISchnittman fc Buonanno! (|2007f ) results were conserva- 
tive: the fraction of large kicks is significantly higher 
than their estimate for comparable-mass mergers with 
plausible spins. 

Barring mechanisms to retain supermassive black holes 
after major mergers, one would expect tens of percent of 
merged galaxies to have no central bla ck hole, in strong 
contr adiction with observations (see iFerrare se fc Ford' 
l2005h . Low spin magnitudes would lower kicks, but 
this is contrary to spin inferences f r om Fe Kg lines; see 
llwasa.wal (119961) :lFabian et al.l(l2002l):lRevnolds fc NowakI 



(|2003f) : iBrennemanfc Reynolds! (|2006D . Ahgnment of 
spins is another possibility; since pure gravity does not 
do this (jSchnittman. ,2004 ; .Bogdanovic et al.. ,2007,1 , ex- 



ternal to rques such as those from nuclear gas would be 
required (jBogdanovic et al]|2007f ). 

In conclusion, we have performed a systematic study 
of the mass ratio dependence of the out-of-plane kicks 
produced by the merg er of spinning black h oles. Our 
work shows that the .Campanelli et ahl (j2007( ) candidate 
kick formula overestimates the out-of-plane kick system- 
atically. However, we find that an additional factor of 
Arj agrees with our numerical results to within 10% (and 
typically 1%) for mass ratios between 1 and 1/3. This 
has considerable implications for black hole retention in 
early dark matter halos, galaxies, and globular clusters. 
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TABLE 3 

Fraction of kick speeds a bove a gi ven threshold, comparei) with the results of 

IScHNiTTMAN &: BuoNANNOl (120071) (SB) AND ICampanelli et al.I (120071) (CLZM). In all cases we 
assume an isotropic distribution of spin orientations. 



Mass ratio and spin Speed threshold SB CLZM This work 

1/10 <<?< 1, ai = a2 = 0.9 > SOOkms^i 0.12(+0.06, -0.05) 0.364±0.0048 0.2283±0.0014 

t;> 1000kms-i 0.027{+0.021, -0.014) 0.127±0.0034 0.085±0.0008 

1/4 < g < 1, ai = 02 = 0.9 i'> 500kms-i 0.31(+0.13, -0.12) 0.699±0.0045 0.618±0.0014 

> lOOOkms-i 0.079{+0.062, -0.042) 0.364±0.0046 0.2547±0.0013 
1/4 < g < 1, < ai,a2 < 1 i'> 500kms-i ■ ■■ 0.428±0.0045 0.3484±0.0015 

> lOOOkms-i ■ ■■ 0.142±0.0034 0.0974±0.0009 



